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Abstract. We study the probability measure /xq for which the moment sequence is 
Cn)n^- ^^ prove that /zq is absolutely continuous, find the density function and 
prove that /Iq is infinitely divisible with respect to the additive free convolution. 



^^ 1. Introduction 

^ A 2-plane tree is a planted plane tree such that each vertex is colored black or white 

. and for each edge at least one of its ends is white. Gu and Prodinger [3] proved, that the 

Ph number of 2-plane trees on n + 1 vertices with black (white) root is ( "^ ) ^^^ (Fuss- 

'^ Catalan number of order 3, sequence A001764 in OEIS |lOj) and (^"^^) 3^^ (sequence 

■^ A006013 in OEIS) respectively (see also [1]). We are going to study the sequence 

^ ^ /3n\ 2 _/3n + l\ 1 /3n + 2\ 2 



m 



X 
S3 



P^ \n J n + 1 \ n J 3n + 1 \ n J 3n + 2 

n, 

which begins with 

t^ 2, 3, 10, 42, 198, 1001, 5304, 29070, 163438, . . . , 

, G^ of total numbers of such trees (A007226 in OEIS). 

Both the sequences on the right hand side of ([I]) are positive definite (see [3, E]), 
^ therefore so is the sequence (^^) ^^ itself. In this paper we are going to study the 

■^ corresponding probability measure /zq, i-e. such that the numbers (^") ^^ are moments 

^ of /io- First we prove that /Xq is Mellin convolution of two beta distributions, in particular 

\Q fiQ is absolutely continuous. Then we find the density function of /iq. In the last 

■^ section we prove, that /io can be decomposed as additive free convolution fii ffl /i2 of 

O two measures, which are both infinitely divisible with respect to ffl and are related to 

the Marchenko-Pastur distribution. In particular, the measure //q itself is ffl-infinitely 

divisible. 



2. The generating function 
Let us consider the generating function 

2^" 






, -^ , n + 1 

n=0 ' 
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According to ([I]), G is a sum of two generating functions. The former is usually denoted 
byiSs: 



n— n \ / 



.„ , 3n + l 

n=0 ^ ' 

and satisfies equation 

(2) Bs{z) = l + z-Bs{zf. 

Lambert's formula (see (5.60) in [2]) implies, that the latter is just square of B^: 






n / 3n + 2' 

n=0 ^ ' 

SO we have 

(3) G{z) = B:i{z) + B^{zf. 
Combining ^ and (|3|, we obtain the following equation for G: 

(4) 2-z-{l + 2z)G{z) + 2zG{zf - z^G{zf = 0, 

which will be applied later on. 

Now we will give formula for G{z). 

Proposition 2.1. For the generating function of the sequence (Til) we have 

(5) G(.) = ""f°;^, . 

(4cos^ a — 1) 

where a = | arcsin f a/272;/4 j . 

Proof. Denoting (a)„ := a(a + 1) ... (a + n — 1) we have 

2(3n)! _-2(f)„+,(^)„+,2r 



n+l 



(n + l)!(2n)! 3(n + 1)! (-i)^_^^ 4- 

Therefore 

9 9 p / -2 -1. 1 I 272^ 

Now we apply formula 



+1 



^- - -1-1 
' 3 ' 3 ' 2 



M I = -i/usin j - arcsin (v^) ) + Vl — m cos I - arcsin (v^) ) 



which can be proved by hypergeometric equation (note that both the functions w ^ 
w sin (I arcsin (w)), w i-)- cos (| arcsin (w)) are even, so the right hand side is well 
defined for \u\ < 1). Putting a = ^ arcsin (a/u), u = 27z/A, we have i/m = sin 3a, 
yj\ — u = cos 3a, which after elementary calculations gives ^5\. D 
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3. The measure 

In this part we are going to study the (unique) measure /iq for which | (^^) ^^ } _ is 
the moment sequence. We will show that fiQ can be expressed as the Mellin convolution 
of two beta distributions. Then we will provide explicit formula for the density function 
V{x) of fiQ. 

Recall (see |llj), that for a,f3 > 0, the beta distribution Beta(a, /3) is the absolutely 
continuous probability measure defined by the density function 

^"•^^"^^ r(a)r(/3) ^ i^ ^^ ' 

for X G (0, 1). The moments of Beta(Q;, /3) are 



r(a)r(a + f3 + n) \_^ a + (3 + i 

For probability measures ui, V2 on the positive half- line [0, oo) the Mellin convolution 
is defined by 



(6) (z/i o z/2) (A) := / / XAixy)diyi{x)du2iy) 

Jo Jo 

for every Borel set A C [0, 00) {xa denotes the indicator function of the set A). This is 

the distribution of the product Xi ■X2 of two independent nonnegative random variables 

with Xj ~ z/j. In particular, if c > then z/ o 5^ is the dilation of u: 

(z/ o 4) (A) = D,z/(A) := z/ (-^A 

where 6c denotes the Dirac delta measure at c. 
If both the measures ui, V2 have all moments 



Sn{j^i) ■= / X"" dl'i{x) 

Jo 
finite then so has Vi o z/2 and 

Sn {Vl O Z/2) = S„(z/i) ■ S„(z/2) 

for all n. The method of Mellin convolution has been recently applied to a number of 
related problems, see for example |6], |8]. 

Now we can describe the probability measure corresponding to the sequence ( ^) ^^. 

Proposition 3.1. Define /xq (^-s the Mellin convolution 

(7) /io = Beta(l/3, 1/6) o Beta(2/3, 4/3) o ^27/4. 

Then the numbers (^)^^ o,re moments of ixq: 



Proof. It is sufficient to check that 



n — 1 ^ i„ , . n—1 



(3n)! _tt1/3 + ^ fr 2/3 + ^ /27 



(n + l)!(2n)! i^^i/2 + ^ IJ 2 + z V4 



AAl/2 + z Ai 

D 
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In view of formula ([T]), the measure fiQ is absolutely continuous and its support is the 
interval [0, 27/4]. Now we are going to find the density function V{x) of yUo- 



Theorem 3.2. Let 

V{x) = 

+ 



^ 



2io/%a;2/3 
1 



h,^l- Ax/27 - l) (l + v^l - Ax/21 
{?,^l-Ax/21 + l) (l + Vl - Ax/27 



28/37^x1/3^3 

X G (0,27/4). Then V is the density function of fiQ, i.e. 



1/3 



-1/3 



27/4 



x" V{x) dx 



3n 



1 



n J n + 1 
for n = 0,1,2, .. .. 

The density V{x) of /io is represented in Fig. l.B. 
Proof. Putting n = s — 1 and applying the Gauss-Legendre multiplication formula 

m — 1 



T(mz) = (2nY^-"'^^''m"''~'^^T(z)T i z + -] T ( z + - 

m I \ m 



T\z + 



m 



we obtain 



3n 



r(3n + l) 



r(3s-2) 



nJn + 1 T{n + 2)T{2n + l) r(s + l)r(2s - 1) 
2 /3"/^27Vr(s-2/3)r(s-l/3) 

~ 27V^ 



Hs) 



A J r(s-i/2)r(s + i) 

Then ip can be extended to an analytic function on the complex plane, except the points 
1/3 -n, 2/3 -n, n = 0,1,2,.... 

Now we are going to apply a particular type of the Meijer G-function, see [9] for 
details. Let V denote the inverse Mellin transform of tp. Then we have 

V{x) = — / x~''ip{s)ds 

27T1 Jc-ioc 

_ 2 fs 1 
~ Tl\n2^i 

2 



c+ioo 



r(s-2/3)r(s-l/3) [Ax\ 



'3 ^2,0 /4^ 
27 V TT ^^'^ I 27 



r(s-i/2)r(s + i) \2i J 

-1/2, 1 



ds 



-2/3,-1/3^ 

where x G (0, 27/4) (consult [11] for the role of c in the integrals). On the other hand, 
for the parameters of the G-function we have 

(-2/3 - 1/3) - (-1/2 + 1) = -3/2 < 

and hence the assumptions of formula 2.24.2.1 in [H] are satisfied. Therefore we can 
apply the Mellin transform on V{x): 



27/4 



x" V{x) dx 



27 V TT 



27/4 



x'~'G. 



1^2,0 
2.2 



Ax 

27 



-1/2, 1 
-2/3,-1/3 



dx 



2 /3 /27 

27V ^ It 



u'~'G 



lr^2,0 
2,2 



u 



-1/2, 1 
-2/3,-1/3 



du = ipi^s) 
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whenever 3fJs > 2/3. Consequently, V = V. 

Now we use Slater's formula (see [9], formula 8.2.2.3) and express V in terms of the 
hypergeometric functions: 

-2/3 



+ : 



27 V vr r(l/6)r(5/6) V27 



r(-l/3) /4x\ 



-2 5 2 



3 '6'3 



Ax 

27 



27 V vr r(-l/6)r(4/3) \27 J 



-1/3 



2i^l 



-17 4 
T'6'3 



4x\ 



x/3 

47rx^/^ 



iFi 



-2 5 2 
ir'6'3 



4x 

27 



+ 



27rV3xV3 



2i^i 



27 J 
-17 4 
3~'6'3 



4x\ 



27; 



Applying the formula 

^ , t-2 t+1 

2-ri 



; t 



2t 



(t - 1 + v^r^^) (1 + VT^~^) 



i-t 



(see [H]) for t = 2/3 and t = 4/3 we conclude the proof. 

4. Relations with free probability 



D 



In this part we are going to describe relations of /io with free probability. In particular 
we will show that /io is infinitely divisible with respect to the additive free convolution. 

Let us briefly describe the additive and multiplicative free convolutions. For details 
we refer to [12], [7] . 

Denote by A^^ the class of probability measures on M with compact support. For 
/i G A^^, with moments 

s„(/i):= [ t^'dfxit), 



and with the moment generating function: 



M,iz) := Y, Uf^)^"" = { 

m=0 ^^ 



dfi{t) 



tz' 



we deflne its R-transform R^{z) by the equation 

(8) R^{zM^{z)) + l = M^{z). 

Then the additive free convolution of fi', fi" G A^^ is deflned as the unique fi'Sfi" G Ai"^ 
which satisfles 

R^i'mf,"{z) = R^'{z) + Rf,"{z). 

If the support of /x G A^'^ is contained in the positive halfline [0, +oo) then we deflne 
its S -transform S^{z) by 



(9) 



M„ 



1 + z 



-SJz] 



l + z 



or Rf,{zS^{z)) 



on a neighborhood of 0. If /i', fi" are such measures then their multiplicative free convo- 
lution fi' Kl /i" is deflned by 

S^'Sfl"{z) = Sf,f{z) ■ S^n{z). 

Recall, that for dilated measure we have: M-d^^{z) = M^{cz), Rji^^{z) = R^{cz) and 
S-D^p{z) = S^{z)/c. The operations ffl and Kl can be regarded as free analogs of the 
classical and Mellin convolution. 
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For t > let Wt denote the Marchenko-Pastur distribution with parameter t: 

r .r jAt-{x-l- ty , 

(10) Wt = maxjl — t, Ojoo H ?; dx, 

znx 

with the absolutely continuous part supported on [(1 — V^)^, (1 + v^)^] • Then 
(U) 

(12) 



^tl^J = 


l + z- 


tz+^l^- 


z-tzf- 


Atz^ 


= 


oo 
n=l 




\k — l) n' 




Rvjt{z) 


tz 


s. 


(z)- ^ 




l-z 


-^"^^ t + 


z' 



(13) 

In free probability the measures Wt play the role of the Poisson distributions. Note that 
from (13) the family {'cut}t>o constitutes a semigroup with respect to ffl, i.e. we have 
lUgSzut = ^s+t for s,t > 0. 

Theorem 4.1. The measure fiQ is equal to the additive free convolution /io = /Ui ffl fi2, 
where Hi = D2U71/2, so that 

/..N 1. V8 - (x - 3)2 



and H2 = nSo + o^i; ^-G- 



1, V4x" 



X 



2 



(15) /i2 = -So H X(OA)i^) dx. 

The measures fj^i, fi2 0,1"^ infinitely divisible with respect to the additive free convolution 
ffl, and consequently, so is ^q. 

The absolutely continuous parts of the measures /ii,/U2 are represented in Fig. l.A. 

Proof. The moment generating function of /iq is M^^^{z) = G{z)/2. Then we have 
M^„(0) = landby Q 

2-z-2{l + 2z)M^,{z) + %zM^,{zf - ^z'M^.izf = 0. 

Let T{z) be the inverse function for Mf^g{z) — 1, so that T(0) = and M^^ (^(-2)) = 
l + z. Then 

2 - T{z) + (-1 - 2T(z))2(l + z) + 8T(z)(l + 2)^ - 8T{z)\l + 2)=^ = 0, 

which gives 

8(1 + zfT{zy - (8^2 + 12;^ + 3)T{z) + 2z = 
and finally 

_ 8z^ + 12^ + 3 - V9 + 8z 4z 

16(l + 2;)3 ~ 8^2 + 12^ + 3 + V9 + 8^' 

Therefore we can find the S'-transform of yUo- 

l + z _ 8^2 + 12^ + 3 - yoTS^ _ 4(1 + ^) 

^^°^'^- z ^''^~ 16^(1 + ^)2 "8^2 + 12^ + 3 + ^9^8^ 
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and from (|9| we get the i?-transforni: 



^'^^^ 2{l-2z) ■ 

Now we observe that R^o{z) can be decomposed as follows: 



Comparing with (13) we observe that Ri{z) is the /^-transform of Hi = Y)2'^\/2i which 
implies that /ii is ffl-infinitely divisible. 
Consider the Taylor expansion of R2{z): 



n=l ^ ^ n=0 ^ ^ 

Since the numbers ( ^) are moments of the arcsine distribution 



n—S^n 



--X{OA){x)dx, 



the coefficients at the last sum constitute a positive definite sequence. So R2{^) is -R- 
transform of a probability measure /i2, which is ffl-infinitely divisible (see Theorem 13.16 
in [7]). Now using ([s]) we obtain 



,, ,, l + 2;z-v/r^4^ 1 l-v^r^4^ 1 1 

My^Az) = = - H = - H , 

"^'^ ' 4:z 2 Az 2 1 + ^/1^4^ 



Comparing with (11 ) for t = 1 we see that /i2 = ^Sq + ^zui. D 

Let us now consider the measures /ii,/i2 separately. For yUi = 02^1/2 the moment 
generating function is 



^■'-'-i..wi-a...-'^g-"S(:)G-i)^ 



n—k 

5 



n=l A;=l 

SO the moments are 

1, 1, 3, 11, 45, 197, 903, 4279, 20793, 103049, 518859, .... 

This is the A001003 sequence in OEIS (little Schroeder numbers), s„(/ii) is the number 
of ways to insert parentheses in product of n + 1 symbols. There is no restriction on 
the number of pairs of parentheses. The number of objects inside a pair of parentheses 
must be at least 2. 

On the subject of fi2, applying ^ we can find the ^-transform: 

2{l + z) _ 1 + z 1 

^'^ ' {l + 2zf l/2 + z'l + 2z 

One can check, that Y]2+i i^ the S-transform oi \5q + \5i, which yields 

(16) ^^2= Q(5o + ^'5ijK/ii. 

We would like to thank G. Aubrun, C. Banderier, K. Gorska and H. Prodinger for 
fruitful interactions. 



WOJCIECH MLOTKOWSKI, KAROL A. PENSON 




(a) The densities of ^i, fi2 (b) The density of fiQ = fiiS ^2 

Figure 1. The densities of /ii, /i2 and /xq = A'-i ffl /^2 
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